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In this paper we present an experimental study of the long surface wave instability that 
can develop when a granular material flows down a rough inclined plane. The threshold 
and the dispersion relation of the instability are precisely measured by imposing a con- 
trolled perturbation at the entrance of the flow and measuring its evolution along the 
slope. The results are compared with the prediction of a linear stability analysis con- 
ducted in the framework of the depth-averaged or Saint- Venant equations. We show that 
when the friction law proposed in Pouliquen (1999a) is introduced in the Saint- Venant 
equations, the theory is able to predict quantitatively the stability threshold and the 
phase velocity of the waves but fails in predicting the observed cutoff frequency. The 
instability is shown to be of the same nature as the long wave instability observed in 
classical fluids but with characteristics that can dramatically differ due to the specificity 
of the granular rheology. 



1. Introduction 

Natural gravity flows such as mud flows or debris flows can be destructive. Very often 
the material does not flow continuously along the slope of mountains, but develops in a 
succession of surges. These surges are of large amplitude and can be destructive as they 
are able to carry large debris. The breaking of an initially continuous flow to a succession 
of waves is often attributed to an inertial instability that can develop in thin liquid flows 
down slopes. However, while this instability has been extensively studied for the case 
of Newtonian fluids, the formation of long surface waves in particulate flows has been 
much less investigated. The goal of this paper is to clarify the dynamics of the long- wave 
instability for a cohesionless granular material flowing down a rough inclined plane. 

Long-wave free surface instability in gravity flows is a phenomenon common to many 
fluids. For a viscous fluid in the laminar regime, the instability is often called "Kapitza 
instability" after the pioneering work by Kapitza & Kapitza (1949). They have observed 
that a thin film of water flowing along a vertical wall does not remain uniform but deforms 
in a succession of transverse waves propagating down the wall. Since Kaptiza's work, 
many experimental and theoretical studies have been devoted to thin viscous flows down 
inclined planes (see, for instance, the review of |Chang 1994] or lOron" et al. 1997j) . Free 
surface waves are also observed in turbulent flows down slopes, where they evolve into 
a series of bores more or less periodic called "roll waves" HCornish 19341 IDressler 19491 
INeedham &: Merkin 19841 |Kranenburg 1992] ) . These roll waves are of great importance 
in hydraulic for flows in open channels. Similar surges are also observed in flows using 
more complex fluids, e.g. mud flows, gravity currents, flows of particles suspension (see 



2 Y. Forterre and O. Pouliquen 

ISimpson 1997| ). However, few precise experimental studies can be found for these non 
Newtonian fluids. 

The spontaneous formation of free surface waves in fluids having very different Theo- 
logical properties comes from the fact that the instability mechanism does not depend on 
the precise fluid characteristics (see the clear discussion by Smith (1990) in the context of 
viscous film flows). When a fluid layer flows down an inclined plane, a small perturbation 
of the free surface propagates at first order with a phase velocity which is different than 
the mean fluid velocity. If the fluid could adjust instantaneously its velocity to the local 
thickness variation, the wave would propagate without amplification (kinematic waves, 
I Wit ham 1974|l . However, because of inertia, the fluid does not immediately adjust its 
velocity when the wave arrives. This delay can give rise to a positive mass flux under the 
wave leading to the growth of the perturbation. 

From a theoretical point of view, many studies concern the case of laminar New- 
tonian flows. In this case the linear stability analysis of a film flow can be derived 
based on the Orr-Sommerfeld equation (Be njamin T957| lYih 19631 . Am plitude equa- 
tions and precise non linear models have been also recently developed ( |Chang 1 994 
|Ruyer-Quil fc Manneville 2000| > . For non Newtonian complex fluids, exact three dimen- 
sional analysis are often not possible because of the lack of sufficient knowledge of the 
constitutive equations. For those systems, a classical approach for studying the stability 
of thin flows is the shallow water description. Mass and momentum balances are written 
in a depth averaged form (Saint- Venant equations. ISaint-Venant 1871(1 . In this frame- 
work, the rheological characteristics of the fluid are mainly taken into account into the 
expression of the shear stress between the flowing layer and the surface of the plane. The 
shear stress is a viscous force for laminar Newtonian flows (Shkado v 1967|) . a turbulent 
friction given by a Chezy formula for turbulent flows ( |Kranenburg 1992] ), a Bingham 
stress for mud flows l|Liu fc Mei 19 94). It is then possible to study the linear stability 
analysis of the flow and to predict the threshold and growth of the waves. 

Although the free surface instability mechanism is the same for different fluids, the 
precise characteristics of the wave development (threshold, growth rate,..) dramatically 
depend on the rheological properties of the material. A precise experimental study of the 
instability could then serve as a test for the rheological laws proposed for complex fluids. 
This is the spirit of this paper on the wave formation in cohcsionlcss granular materi- 
als. Formation of long waves in granular flows down inclined planes has been reported 
in previous studies USavage 1979| IDavies 19901 IValTa ncc 1994j |Ancey 1997|lDae7r 2001jl . 
However, to our knowledge no precise measurements of the instability properties have 
been carried out. 

The rheology of dry granular materials flowing in a dense regime is still an open problem 
( |Rajchenbach 2000| |Pouliquen fc Chevoir 200"2"1 ). However, for the flow of thin granular 
layers on inclined planes, a Saint- Venant description seems to be relevant. Savage and 
Hutter first proposed this approach for describing the motion of a granular mass down 
an inclined plane {Savage Hutter 1989[). For the basal friction stress they used a solid 
friction law where the shear stress at the interface between the flowing material and the 
inclined plane is proportional to the normal stress, i.e. to the weight of the material col- 
umn above the base. Although this description is able to correctly predict the spreading 
of a mass on steep slopes and smooth surfaces ( Savag ed Hutter 19 89 Naai m et al. 19971 
|Gray et al. 1 999 Wiclan cTet al. 1999}) . it does not capture the existence of steady uni- 
form flows when the surface is rough and the material sheared in the bulk. Recently, 
a generalization of the basal friction law has been proposed based on scaling properties 
observed in experiments ( |Pouliquen 1999a| ) . The generalization gives a friction coefficient 
which is no longer constant but depends on the local velocity and thickness. Using this 
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law in Saint- Venant description it is possible to quantitatively predict the motion of a 
mass down a rough inclined plane from initiation to deposit ( |Pouliquen fc Forterre 2002"| ). 
In this paper we want to check to what extent the same approach using the generalized 
friction law can quantitatively predict the long wave instability for dry granular materials. 

To this end we have carried out precise measurement of the linear development of 
the instability in the same spirit as Liu et al. (1993) did for liquid films. Liu et al. 
have developed a method to induce at the entrance of the liquid flow a perturbation 
whose amplitude and frequency can be controlled. The instability being convective, the 
observed waves downstream result from the amplification of the entrance perturbations. 
By following the evolution of the injected perturbation for different frequencies, Liu et 
al. were able to experimentally measure the dispersion relation and precisely determine 
the instability threshold. In this paper we follow the same procedure for granular flows: 
in the limited range of inclinations where steady uniform flows are observed, we impose 
an external forcing to control the wave development. 

The paper is organized as follows. We first present preliminary experimental observa- 
tions of the instability without forcing (§2). We show that a simple visual analysis of 
the free surface instability is not sufficient to conclude about the instability mechanism 
and that more precise measurements are necessary. The experimental set up, the forcing 
method and the measurement procedure are described in § 3. The linear stability analysis 
based on Saint- Venant equations is presented in § 4. Results and comparison between ex- 
periments and theory are presented for glass beads in § 5 and for sand in § 6. Discussion 
and conclusion are given in § 7. 

2. Preliminary observations 

Free surface waves in granular flows down a rough inclined plane have been reported 
previously by several authors. Using glass beads, Savage (1979) and Vallance (1994) ob- 
served the deformation of the free surface in a succession of surges down a long chute flow. 
Ancey (1997) reported formation of surge waves with sand. Recently, Daerr (2001) has 
noticed an instability at the rear of avalanche fronts for flows of glass beads on velvet. All 
these experiments were performed on rough surfaces in a dense regime, i.e. with a mean 
volume fraction of the medium almost constant. Waves have also been recently observed 
for granular flows on flat smooth surfaces l|Prasad et al. 2000"l|Louge fc Keast 2001| | but 
dramatic variations of density are observed in this case, with very dilute regions. We will 
not discuss this case here. 

FigureHshows the typical free surface we observe for a flow of sand 0.8 mm in diameter 
down an inclined plane made rough by gluing one layer of sand grains on it. A silo at the 
top of the plane (2 m long, 70 cm wide) continuously supplies the flow at a constant flow 
rate. Waves first appear close to the outlet of the silo as two-dimensional deformations. 
They rapidly amplify downstream and break down in the transverse direction due to 
secondary instabilities. The shape of the saturated waves is highly non linear as shown 
by the typical thickness profile plotted in inset of figure ^ The distance between two 
surges is typically of the order of 20 cm, which is very large compared to the thickness 
of the layer (which is about 5 mm for the case of figure QJ. 

The observed pattern does look similar to patterns observed in liquid flows. However, 
a major difference appears when changing the experimental conditions, i.e. the thickness 
of the layer and the inclination of the plane. Figure [2] shows a qualitative phase diagram 
obtained simply by visual observations of the free surface. The striking result of this pre- 
liminary analysis with sand is that the waves are observed mainly close to the threshold 
of the flow (low inclinations and thin flows). In this region, large amplitude waves are 
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Figure 1. Free surface of the sand flow showing the long surface waves. Material is flowing 
downwards; the plane width is 7 cm. Inset: typica l thickness profile of the natural waves mea- 
sured with a laser sheet light (see |Pouliquen 1999a| for the description of the method). The arrow 
gives the flow direction, (sand 0.8 mm in diameter, 6 = 34°, h =4.6 mm). 

observed. When one increases the thickness of the layer at a fixed inclination, the waves 
form further and further downstream and eventually disappear for thick flows i.e. for 
more rapid flows. Naively, one would say that this observation is in contradiction with 
the explanation in terms of an inertial instability, which should develop preferentially in 
rapid flows rather than in slow flows. 

Another striking preliminary observation is made using glass beads as a granular mate- 
rial. In this case, no wave was observed in our set-up, whatever the thickness and the incli- 
nation. The absence of noticeable deformation of the free surface is what made systematic 
measurements of the steady uniform regime possible and allowed to put in evidence flow 
scaling laws j |Pouliquen 1999a| ). This point is also intriguing as one would expect the 
long- wave instability to occur as in other fluids for fast enough flows (JVardoulakis 2002 ) . 

The situation is then rather confused and the preliminary observations raise several 
questions. First, why are the flow of sand and the flow of glass beads so different? Sec- 
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Figure 2. Qualitative stability diagram for sand 0.8 mm in diameter. The hatched zone is the 
domain where waves are observed. (•): waves, (o): no wave, (□): flow threshold. 

ondly, in sand the waves seem to develop for slow flows and disappear for rapid flow 
in apparent contradiction with an inertial instability. Are the waves due to a different 
instability mechanism? Finally, why is no wave observed with glass beads whereas for 
rapid flows one would expect the instability to develop? In order to properly answer those 
questions one need to precisely study the stability of the flow not only by visual obser- 
vation of the free surface deformation but by careful measurement of the amplification 
of an initially imposed perturbation. This is the reason why, following Liu and Gollub's 
procedure (Liu et al. 1993), we have developed a forcing method. 

3. Experimental setup 

The experimental set-up is presented in figure E] The plane is 2 m long and 35 cm 
wide. The bottom plate as well as the side walls are glass plates. The roughness is made 
by gluing on the bottom plate one layer of the particle used for the flow. The flow rate is 
controlled by the opening of the silo. In this paper we have used two granular materials: 
glass beads d — 0.5 mm in mean diameter (the same as used in |Pouliquen 1999a| ), and 
sand d — 0.8 mm in mean diameter (figure^). 

3.1. Forcing method 

In order to impose a perturbation whose frequency and amplitude can be easily con- 
trolled we periodically blow a thin air jet on the free surface. The jet is created by three 
loudspeakers embedded into a two dimensional nozzle with a 1 mm slit (figure Et). A 
homogeneous and localized air jet is then created, with an amplitude and a frequency 
controlled by the amplitude Aqbf and frequency / sent to the loudspeakers via a low 
frequency function generator. The nozzle is placed 30 cm from the outlet of the reservoir. 
The typical free surface deformation we achieved with this set up is of the order of 0.25 
mm and frequency varies between 1 Hz and 20 Hz. A typical response of the free surface 
to the forcing air jet is plotted in figure Eta). In this figure, the measurement is carried 
out just below the nozzle and the signal sent to the loudspeaker is sinusoidal at a fre- 
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Figure 3. (a) Sketch of the experimental set-up used to force the instability, (b) Description 

of the light detection method. 



quency of 3 Hz. One observes that the thickness variation follows the forcing frequency 
but that the response is not sinusoidal but highly asymmetric. When air is blown, the 
air jet is localized and the induced deformation important. When the nozzle sucks the 
air, there is no influence on the free surface. This difference between localized ejection 
and diffuse injection is well-known and is used, for instance, for propulsion in water. The 
non sinusoidal response clearly appears on the Fourier transform of the signal in figure 
[5^6). The perturbation is periodic at the imposed frequency, but many harmonics are 
also present with amplitude that can be higher than the fundamental. This effect is more 
pronounced when the forcing frequency is low. With this forcing method it is then diffi- 
cult to inject low frequency modes at measurable amplitude, without injecting very high 
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Figure 4. Granular materials used, (a) Glass beads d = 0.5 ± 0.05 mm in mean diameter 
(grinding glass beads supplied by Potters-Ballotini, France), (b) Sand d = 0.8 ±0.1 mm in 
diameter (sand of the Loire supplied by Sifraco, France). 
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Figure 5. Free surface deformations of the flow below the nozzle (/ =3 Hz, Aobf = 1-5 V). 
(a) Thickness variations, (b) Fourier transform of the thickness variations, |a/|, averaged over 
60 cycles (arbitrary units). 

harmonics, which will interact in a non linear regime. Consequently, no measurement of 
the linear evolution was possible below 1 Hz. 



3.2. Thickness measurement 

In order to measure precisely the free surface deformation we use a light absorption 
method. A parallel beam lights the plane from above as sketched in figure Two 
photodetectors are placed at different distances from the nozzle below the plane. The 
light measured by the detector varies exponentially with the thickness of the granular 
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Figure 6. Calibration curve for the light transmitted through a layer of grain as a function 
of the thickness h of the layer, (a) glass beads d = 0.5 mm in diameter, (b) sand d = 0.8 mm 
in diameter. (•): static layers, (o): flowing layers. The solid curve is the best exponential fit 
I/Io = exp(—h/L a ) with L a = 7.40 d for glass beads and L a = 1.82 d for sand. 



layer as shown by the calibration curves in figure EJ The light intensity received by the 
detector is plotted versus the thickness (see |Pouliquen 1999o| for the thickness measure- 
ment method). The data are well fitted by I = Iq exp(— h/L a ), the attenuation length 
being larger for the beads L a —7A d, than for the sand i a =1.82 d. This calibration has 
been carried out for both uniform layer at rest (black dots) or flowing layers (circles). 
We observe that the two data sets collapse. Since the light absorption is a function of 
both the layer thickness and the volume fraction, this collapse means that the variations 
of the volume fraction are negligible in the dense flow regime studied here. 

The signal measured by the photodetector is digitalized by a PC trough an acquisition 
board at 100 kHz. Notice that from the temporal evolution of the light amplitude, J(t), 
it is possible without any calibration to get the deformation of the free surface 8h(t) = 
h(t) — (h(t)) where (h(t)) is the mean thickness. The exponential attenuation law leads to 
the following expression: 8h(t) = —L a (\nl(t) — (\nl(t))), which is independent of Jo and 
of the gain of the photodetectors. With this method the deformation of the free surface 
is determined with a precision of about 0.05 mm. 



3.3. Measurement of the dispersion relation 

From the thickness measurement at two different points it is in principle possible to 
measure the dispersion relation. To this end a perturbation at a given frequency / is 
imposed to the flow by our forcing device. The two photodetectors then give the free 
surface deformations 5h{x\,t) and 5h{x2,t) at two different locations x\ and X2- By 
computing the Fourier transform of the deformations, one obtains at the two locations 
the amplitudes cif(xi) and 0/(2:2) and the phases (f>f(xi) and 4>f i x 2) of the fundamental 
mode at frequency / defined by : Shf(x,t) = a f (x) cos(2ir ft — <fif(x)). If the evolution 
between x\ and X2 is assumed to be exponential, i.e. in the linear regime of the instability, 
the deformation of the mode / is supposed to vary as: 5hf — aoe ax cos(2nf(t — x/v^)). 
The spatial growth rate a(f) and the phase velocity v,p(f) are then given by: 



Long surface wave instability in dense granular flows 



9 



-(/) = ^rln('?4^V (3-D 



x 2 - x\ \df(xi) 
2n{x 2 -x l )f 

Mf) = ~T~t — \ J 1 \ ' ( 3 - 2 ) 

The dispersion relation can then in principle be determined. Experimentally the diffi- 
culty is to determine a region of linear instability where the wave grows exponentially. 
The amplitude of the deformation has to be small enough to remain in the linear evolution 
but not too small to be measured. 



4. Linear stability analysis 

In this section, we present the theoretical analysis of the long-wave instability for 
granular flows in the framework of the depth-averaged equations. A previous analysis has 
been carried out by Savage (1989) using a Chezy formula to model the basal shear stress. 
Here we present the theoretical prediction when the interaction between the granular layer 
and the rough plane is described by the empirical friction law derived from experimental 
measurements on steady uniform flows jPouliquen 1999a). 

4.1. Depth-averaged equations 

A detailed derivation of depth-averaged equations in the context of granular flows can 
be found in Savage & Hutter (1989). Assuming that the flow is incompressible and that 
the spatial variation of the flow takes place on a large scale compared to the thickness 
of the flow, we obtain the depth-averaged mass and momentum equations by integrating 
the full three-dimensional equations. For a two-dimensional flow down a slope making an 
angle 8 with the horizontal (see figure 0), depth-averaged equations reduce to: 



dh dhu 

()l ()r j = (ta,n6-iJ,(u,h)-K^)pghcos6, (4.2) 

where h is the local thickness of the flow and u is the averaged velocity defined by 
u = Q/h, Q being the flow rate per unit of width. 

The first equation l|4.1|) is the mass conservation. The second equation (|4.2|) is the mo- 
mentum equation where the acceleration is balanced by three forces. In the acceleration 
term, the coefficient a is related to the assumed velocity profile across the layer and is of 
order 1. We will discuss the role of the parameter a in § 5. The first force in the right-hand 
side is the gravity parallel to the plane. The second term is the tangential stress between 
the fixed bottom and the flowing layer; it is written as a friction coefficient /i multiplied 
by the vertical stress pghcosO. The friction coefficient is assumed to depend only on the 
local thickness h and the local velocity u. This is a generalization of Savage & Hutter 
(1989) 's assumption of a constant solid friction. However, the friction coefficient p could 
depend on other quantities such as the normal stress or the derivatives of h and u. The 
last term in (|4.2|) is a pressure force related to the thickness gradient. The coefficient K 
represents the ratio of the normal stress in the horizontal direction (cc-direction) to the 
normal stress in the vertical direction (z-direction). Recent numerical simulations have 
shown that for dense granular flows the horizontal and the vertical normal stresses are 
almost equal (Prochnow et al. 2000, Ertas et al. 2001). Therefore, we will take K = 1 



10 



Y. Forterre and O. Pouliquen 




in the following. However, it should be noticed that for non uniform flows, the factor K 
could depend on the divergence of the flow QGray et al. 1999| IWieland et al. 1999|l . 

The great advantage of the Saint- Venant equations is that the dynamics of the flow- 
ing layer can be predicted without knowing in details the internal structure of the flow, 
although some information of the flow dynamics is lost. The complex and unknown 
three-dimensional rheology of the material is taken into account only in the friction term 
/i(u,h). Moreover, this friction is easily determined experimentally by studying steady 
uniform flows. These flows simply result from the balance between gravity and friction, 
i.e. fx(u,h) = tan#. This means that measuring how the mean velocity u varies with 
the thickness h and inclination 9 is sufficient to determine the friction law. Knowing the 
function u(h, 6) is then equivalent to knowing the friction law (see |Pouliquen 1999a). 





4.2. Stability analysis 

We study here the stability of a steady uniform flow of thickness h and averaged velocity 
Uq. For the sake of simplicity, we first re- write the Saint- Venant equations (|4.1|l and 
(|4.2I) with the dimensionless variables given by: h — h/h , u ~ u/uq, x — x/ho and 
t = (uo/ho)t. We then obtain: 

(4.3) 

= ( t&n6-Li(u,h) - ^ ] h, (4.4) 
where F is the Froude number defined by: 

yglio cost/ 

The two dimensionless control parameters of the problem arc therefore the Froude number 

F and the angle of inclination 8. 

We look for steady uniform solutions given by: 

h(x,t) = 1, u(x,t) = 1. (4.6) 

For this basic state, the mass conservation <|4.3ll is satisfied and the momentum equation 
(|4.4|) reduces to the balance between gravity and friction: 1) = tan (9. The next 
step is to study the stability of 1)4.6(1 by perturbing the flow: h(x,t) = 1 + hi(x,t), 
u(x,t) = 1 + u\(x,t) with (hi(x,t) <§; 1, u\(x,t) <C 1) and by linearizing the depth- 
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averaged equations which become: 

§ + ^ + ^=0, (4.7) 
at ox ox 

(it + {a - + {2a ~ 1) ii) = - 6/11 - 1£> (48) 

where the dimcnsionlcss variables a and & are related to the friction law fi(u, h) by: 



(the index '0' means that the derivatives are calculated for the basic state). Note that to 
derive l|4.8(l . we have used the mass conservation 14.7(1 . 

We then seek normal mode solutions for the perturbations: hi(x,t) = Ziexp i(kx — Cot) 
and ui(x 7 t) = uexpi(fci — Cot) where k is the dimensionless wavenumber and Co is the 
dimcnsionlcss pulsation, which, once introduced in the linearized equations l|4.7f) and 
(|4.8II . gives the following dispersion relation: 

_ 2 r,~r * /- . * ~ -A / 1 
— w + zauik 



F 2 



((a-b)k-ati) + (j^-c^Jk 2 = 0. (4.10) 



From this dispersion relation, we can easily compute the spatial growth rate and the 
phase velocity of the waves for a given imposed real pulsation. Computation of the 
spatial stability analysis is given in Appendix A. This analysis shows that the flow is 
unstable when: 



l-->a+\a(a-l) + -^. (4.11) 
a V r 

The above stability criterion is expressed as a function of a and b written in terms of 
the derivatives of the friction law fj,(u, h). Experimentally, we have access to the relation 
u(h, 9) for steady uniform flows. It is therefore useful to write a stability criterion using 
the relation u(h, 9) instead of h). Moreover, we will see that this formulation gives 
an interpretation of the long-wave instability in terms of waves interactions. 

As for steady uniform flows we have n(u, h) = tan(#(it, h)), it is easy to show that: 

1 (d6\ b fdu\ 

a -c^N ' and a=-UJo' (412) 
The stability criterion l|4.11|l may therefore be written as: 

co > c+, (4.13) 

where 

* = !+(%), (4-14) 



dhj o 



5+ =a+^ja(a-l) + ^. (4.15) 

It can be shown that cq is the (dimensionless) velocity of the kinematic waves, whereas 
c+ is the (dimensionless) velocity of the "gravity waves" that propagate downstream (see 
Appendix B). The stability criterion l|4.13|l therefore means that the flow is unstable 
when the velocity of the kinematic waves is larger than the velocity of the gravity waves. 
Using this formulation, one can then conclude about the stability of the flow just from 
the velocity law u(h, 9) of steady uniform flows. 
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Figure 8. (a) Measurements of the averaged velocity of steady uniform flows u as a function of 
the thickness h and the angle of inclination 9 for glass beads d =0.5 mm (o) UPouliquen 1 999a I 
and sand d =0.8 mm (•). The data collapse on a curve when u/yfgK is plotted as a function of 
h/h s top{0). (b) deposit function h 3 t op (8) for the glass beads (o) and the sand (•). The solid line is 
an interpolation given by h stop (9) = L((tanS2 — tan<5i)/(tan 6 — tan 5\) — 1) where (L/d = 1.65, 
Si = 20.90°, S 2 = 32.76°) for the glass beads and (L/d = 2.03, Si = 27.0°, S 2 = 43.4°) for the 
sand. 



4.3. Velocity law for granular materials 

In this study we use two different granular materials: glass beads and sand. The velocity 
law for the glass beads has been previously measured l |Pouliquen !!)!)!)(/;. The mean 
velocity u, the inclination 8 and the thickness h are related through the following relation: 

P~L — ~7a\- ( 4 -!6) 



\fgh h stop (9) 

where (3 — 0.136. The function h s top{0) is the thickness of the deposit left by a steady 
uniform flow at the inclination 9 jPouliquen fc Renaut 1996||Daerr fc Douady 1999| ). 

In order to know the velocity law for the 0.8 mm sand we have also performed the same 
systematic measurements of the mean velocity as a function of thickness and inclination. 
We observe that, as for glass beads, the velocity is correlated to the deposit function 
h stov (0) as shown in figure |SJ However, the master curve is different for sand and glass 
beads. 

We can therefore write the velocity law in the cases of glass beads and sand in a similar 
form: 

-7 + /?--^ (4-17) 



Vgh h st op(0) 
with 

7 = (5 = 0.136 for glass beads, 

7 = 0.77 (3 = 0.65 for sand, 

(the function h stop (6) for the glass beads and the sand is given in the caption of figure 

lEJ. 

From this relation we can then apply the stability criterion (|4.13(l which compares the 
kinematic wave velocity and the gravity wave velocity. To compute the velocities, one 
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needs to know the coefficient a and b given by the (|4.12(l . From Ij4.17|l . one obtains: 

— H ^ 6) F and*=-^ + 7 \ (4.18) 



cos? / 2 6h s top {0) (FVcos0 + 7)' a \2 Fy/cos . 
The (dimensionless) velocity of the kinematic waves given by l|4.14|l is therefore: 

^ = \ + T7=» (419) 

l Fy cos 

The speed of the gravity waves depends on the parameter a l' see I4.15)) . With a — I, the 
instability condition takes a particularly simple expression: 

The linear stability analysis predicts therefore an instability for granular flows when the 
Froude number is above a critical Froude number F c . For glass beads (7 = 0) the critical 
Froude number is independent of the angle of inclination and given by F c = 2/3. In the 
next section, we will compare these theoretical predictions to the experimental studies. 



5. Results for glass beads 

We first present the experimental results for the flow of glass beads. When there is no 
forcing, no wave was observed. By imposing a forcing at the entrance of the flow we show 
that the instability exists but that the spatial growth rates are small. Our inclined plane 
was therefore too short to allow the observation of the natural instability. 

5.1. Spatial evolution of a forcing wave: existence of a linear regime 

An example of an amplification of the imposed perturbation is shown in figure OH The 
typical temporal signals measured at two locations in the unstable regime (6 = 29°, h = 
5.3 mm) is shown for a forcing frequency / —2 Hz. First, we note that the instability is 
convective, i.e. the periodic perturbation imposed at the entrance is carried downstream 
by the flow. We also observe that the wave strongly evolves along the plane: in the power 
spectra figuresEH&) and0(d) the fundamental mode and the first harmonic of the forced 
wave are amplified downstream whereas high-frequency harmonic modes are damped. 
The amplification of the low-frequency perturbations is also observed on the spectra of 
the natural noise. 

At this stage, an important question is whether the amplification of the low-frequency 
modes of the forcing waves results from a linear instability or from non linear interactions 
between modes. FigureUHlshows the spatial evolution of the fundamental mode of a forced 
wave along the slope for different forcing frequencies (the control parameters are the same 
as in figure EI). We observe that the forcing modes evolve exponentially along the plane, 
which implies that we are studying the linear development of the instability. 

Another proof of the existence of a linear region is given in figure [TIT a) . which shows 
the spatial growth rate a of a forced mode as a function of its frequency /. The black 
bullets are measurements obtained by imposing a forcing at different frequencies / and 
measuring the growth rate of the fundamental mode /. The other curves are obtained 
by keeping the forcing frequency /o constant and measuring the growth rate of all the 
harmonics /o, 2/o, 3/o- . . . We observe that both methods give the same curves cr(f), 
which shows that the different harmonic modes of the forced waves weakly interact and 
evolve quasi-independently in a linear regime. 
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Figure 9. Wave measurement at two distances from the nozzle, x = 6 cm and a; = 96 cm, in 
the unstable regime. (/ = 2 Hz, # = 29°, /i = 5.3 mm), (a) and (c) Local thickness variations. 
(6) and (d) Power spectra averaged on 40 cycles. 

5.2. Dispersion relation 

The existence of a linear region for the instability allows the precise measurement of 
the linear dispersion relation of the surface waves. To measure the growth rate, the two 
photodetectors are separated by 60 cm. 

A typical experimental dispersion relation is presented in figure^JJfor the above control 
parameters (0 = 29°, h = 5.3 mm). First, we note in figure lTTT a) that the spatial growth 
rate of a mode varies with the frequency. For low-frequencies, the spatial growth rate is 
positive, i.e. the amplitude of the mode is amplified along the plane as shown in figure 
1101(2 Hz and 5 Hz). The flow is therefore unstable for the considered control parameters. 
There exists a cutoff frequency for which the amplitude of the mode remains constant 
all along the plane (see the mode / —7 Hz in figure ITU)) . For higher frequencies, the 
amplitude decays along the plane and the mode is stable (see the mode / =12 Hz in 
figure ITU)! . We shall see that the existence of a well-defined cutoff frequency f c will allow 
us to determine precisely the stability threshold. The cutoff frequency f c can typically 
be determined with a precision of about 0.5 Hz. 

The phase velocity of the wave as a function of the frequency can also be determined 
as shown in figure HTT M . We observe that the phase velocity is almost independant 
of the frequency (the variations of are within the error bars) , which means that the 
system is non dispersive within the range of frequency explored in the experiment (see 
also the inset). For 9 = 29° and h = 5.3 mm, v$ «45 cm/s, which is a little more than 
the double of the average velocity of the flow (u «22 cm/s). We have systematically 
measured the phase velocity for other values of the control parameters (9,h). The phase 
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Figure 10. Spatial evolution of the amplitude of the fundamental mode for different forcing 
frequencies. The distance x is measured from the nozzle. The solid curves are the best exponential 
fits exp(cra;) with a (2 Hz) =0.0194 cm" 1 , a (5 Hz) =0.0110 cm" 1 , a (7 Hz) =0.0001 cm" 1 and 
a(Y2Hz) =-0.0389 cm" 1 . The error bars indicate the dispersion of the results. (0 = 29°, h = 5.3 
mm). 



velocity is always nearly independent of the frequency. However, its value depends on the 
parameters (0,h) as we shall see in § 5.4. 

5.3. Stability threshold 

So far, the results have been presented for a given inclination 9 and thickness of the flow 
h. When decreasing the thickness of the flow at fixed 6, we observe that both the cutoff 
frequency and the growth rate of the most unstable mode decrease as shown in figure 
1121 The flow is then less unstable as the thickness of the flow decreases and eventually 
becomes stable below a critical thickness h c . 

In order to measure precisely the stability threshold h c , we have systematically mea- 
sured the cutoff frequency f c as a function of the thickness h as shown in figure ITST a'). 
The neutral stability curve f c {h) is then extrapolated to zero frequency in order to de- 
terminate h c . In order to get a good estimation of the critical thickness we have chosen 
to interpolate the stability curve f c (h) with two polynomials (see the caption of figure 
I13|) . For 9 = 29°, we obtain h c (29°) — 2.6±0.3 mm. The same method can be applied to 
measure the critical Froude number F c , where the Froude number, F = u/y/ghcos9, is 
measured for each set of (8,h) (see figure [T3b) . In that case, the two control parameters 
are (9,F) instead of (9,h). For 9 = 29°, we find F c (29°) = 0.54±0.02. 

The entire procedure is applied for different angles of inclination. Results are shown in 
figure^Jwhich gives the cutoff frequency of the waves as a function of the Froude number 
F. Note that no measurement was made for 9 above 29° or 9 below 24°. For 8 > 29°, the 
flow is no longer in a steady uniform regime but accelerates along the slope and leaves 
the dense regime. For 9 < 24°, the flow rate required to reach the instability increases 
dramatically. Even with an amount of granular material as large as 150 liters, the total 
duration of the flow at the threshold is not long enough to allow precise measurements 
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Figure 11. Experimental dispersion relation (6 — 29°, h — 5.3 mm), (a) Spatial growth rate a 
as a function of the frequency measured from the fundamental mode of the forced wave (•) and 
from the harmonic modes (in that case, the forcing frequency is given in the inset). (6) phase 
velocity as a function of the frequency. Inset: pulsation uasa function of the wavenumber k. 



of the stability threshold. This limitation also explains the large uncertainty at 9 = 24° 
and why no measurement has been made for high Froude numbers at this angle. 

We observe in figurc^Jthat the cutoff frequencies obtained at different angles collapse 
close to the threshold when plotted as a function of the Froude number. This result 
implies that the instability is controlled by the Froude number, at least near the stability 
threshold. 

From these measurements we can determine the stability diagram of the flow in the 
phase space (0,h) or (9,F), as presented in figure IT51 We note that the instability takes 
place for high inclinations and high thicknesses. The control parameter is the Froude 
number as shown in figure [rBT 6). While the critical thickness h c strongly varies with the 
angle of inclination of the plane, the critical Froude number F c weakly depends upon 9. 
For glass beads, we find F c — 0.57 ± 0.05. 
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Figure 12. Spatial growth rate as a function of the frequency for different thicknesses of the 
flow h and a fixed angle 6 = 29°. The corresponding Froude numbers are 1.02 (•), 0.66 (■), 0.58 
(o) et 0.51 (□). 




(a) (b) 

Figure 13. Neutral stability curve for 6 = 29°. The cutoff frequency f c is shown either as a 
function of the thickness of the flow h (a) or as a function of the Froude number F (b). The solid 
curve is a polynomial fit of the type F,h = a + bf c + cf% and the dashed curve is a quadratic fit 
of the type F,h = a + cfc ■ 



5.4. Comparison with the theory 

The main experimental results of the instability for glass beads are: 

• The stability threshold is controlled by the Froude number. Waves appear above a 
critical Froude number given by F c = 0.57 ± 0.05. 

• The surface wave instability is a long-wave instability, i.e. the first unstable mode is 
at zero frequency (zero wavenumber). Above the stability threshold, there exists a cutoff 
frequency for the instability. 

• Within the frequency range investigated, the phase velocity of the waves weakly 
depends upon the frequency, i.e. the media is non dispersive. 
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Figure 14. Neutral stability curve f c {F) for different inclination 9. The continuous curves 
shows for each angle a fit of the type F = a + bf^. 



In this section, we compare these results with the prediction of the linear stability 
analysis performed in the framework of the depth-averaged equations (see §4). We have 
seen that this analysis gives the stability threshold and the spatial relation of dispersion 
as a function of the velocity law u(h, 9) which, for glass beads, is given by (|4.1tj|) . 

The only parameter that is unknown in the theory is the parameter a in the acceleration 
term of the Saint- Venant equations. This parameter is related to the unknown velocity 
profile across the layer (see § 4.1). For a granular flow down an inclined plane, there is no 
experimental measurement of the velocity profile but a is close to one. In the following, 
we will take a = 1 for the sake of simplicity. We will discuss in more details the influence 
of a at the end of the section. 

5.4.1. Stability threshold 

We first compare the theoretical and the experimental stability threshold. For glass 
beads, we have seen that the linear stability analysis predicts an instability when: 

F>1> (5.1) 

2 %/ cos 9 , 
or h > - h stop {6). (5.2) 

The theoretical prediction is given in figure ^| by the dashed curve. We note that the 
agreement between the experimental data and the theory is relatively good concerning 
the variation of the critical thickness with the angle and for the order of magnitude of 
the critical Froude number. However, the theory predicts a critical Froude number that 
is about 15 — 20% higher than the experimental one. 
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FlGURE 15. Stability diagram for glass beads in the phase space (9, h) (a) and (9, F) (b). (•): 
experiments, dashed lines: theoretical predictions. The circles (o) in (a) is the deposit function 
hstop' 



5.4.2. Dispersion relation 

The linear stability analysis gives also the spatial dispersion relation for the waves. 
A typical comparison between theory and experiment is given in figure ^3 which shows 
on the same graph the predicted dispersion relation and the experimental data for a 
given set of control parameters (6 = 29°, F =1.02). We note that the main discrepancy 
between theory and experiment is that the linear stability analysis does not predict 
a cutoff frequency for the waves (figure 116b ') . No term in the Saint- Venant equations 
stabilizes the short wavelengths. However, the linear stability analysis gives the good 
order of magnitude for the maximum growth rate of the instability, when one compares 
the maximal growth rate measured in the experiment and the maximal growth rate 
predicted by the theory. 

It is also interesting to compare the experimental phase velocity to the one predicted 
by the linear stability analysis (figure llfib ^l . We observe that within the frequency range 
investigated in the experiment, the theory predicts a constant phase velocity, in good 
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Figure 16. Theoretical (dashed line) and experimental (•) dispersion relation for 8 = 29° and 
F =1.02. (a) spatial growth rate <r as a function the pulsation Cb. (b) phase velocity v^. a, v$ 
and u are nondimensionalized using the thickness of the flow h = 5.3 mm and the averaged 
velocity of the flow u =21.7 cm/s. 




Figure 17. (•) experimental phase velocity as a function of the Froude number F for different 
angles of inclination (24° — 29°). For each (9, F) parameter, is obtained by averaging the 
dispersion relation over the frequency. The solid curve is the velocity of the gravity waves 
c+ = l + (l/F) (q = 1). 

agreement with the experimental data. The theory also predicts an increase of the phase 
velocity for very low frequencies which we are unable to test with our forcing method. 
It must be noticed that the two limits of the theoretical phase velocity have a precise 
physical meaning. For Co = 0, is equal to the velocity of the kinematic waves cq = 5/2 
while for Co — ► oo, tends to the velocity of the gravity waves c+ = 1 + (1/F) (see § 
4.3). 

We have confirmed this result by systematically measuring the phase velocity of the 
waves for different angles of inclination and Froude numbers. The data are presented in 
figure El We observe that the experimental phase velocity is always close to the speed 
of the gravity waves (solid line) . 
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5.4.3. Influence of the parameter a related to the velocity profile 

The results of the linear stability analysis presented previously are obtained with the 
parameter a equal to 1. This parameter appears in the Saint- Venant equations in the 
acceleration (§4.1) and is defined by a = (u 2 )/(u) 2 . It is therefore necessary to make 
an assumption on the velocity profile across the layer to know the value of a. The value 
a = 1, corresponding to a plug flow, was used in the pioneered work of Savage & Hutter 
(1989). More recent studies on granular surface flows assume that the velocity profile is 
linear and use therefore the value a — 4/3 (Khakar et al. 1997, Douady et al. 1999). 
For a granular flow down a rough inclined plane, recent numerical simulations suggest 
that the velocity profile does not remain linear for thick layers but is close to a Bagnold 
profile IjErtas et al 200l|l . i.e. varies as z 3 / 2 (a = 5/4). From the instability criterion 
l|4.13[) . we obtain for a = 5/4 (resp. a = 4/3 ) a critical Froude number F c = 0.89 (resp. 
F c — 1, 03) to be compared with the experimental value F c = 0.57 ± 0.05. Paradoxically, 
the prediction of the theory (for both the threshold and the wave velocity) seems to be 
better with a = 1. 

These results show that taking into account the velocity profile by a constant parameter 
a in the depth-averaged equations is not very satisfying. This result is well known in 
the case of viscous liquid films. It has been shown (see |Ruyer-Quil fc Manneville 2000| 
for instance) that the simple Saint- Venant equations do not predict quantitatively the 
primary instability if one introduces a = 6/5 corresponding to the parabolic velocity 
profile observed in viscous flows. Comparison between the full linear stability analysis 
of the Navier-Stokes equations and the prediction of the Saint- Venant equations has 
shown that the latter overestimates the stability threshold by about 20%. This difference 
between the exact three-dimensional resolution and the averaged equations comes from 
the fact that for non stationary flows, the shear stress at the base is no longer exactly 
given by its expression for steady uniform flows ( |Ruyer-Quil fc Manneville ~0 00 ) . 

It is therefore doubtful to assign a precise physical meaning to the value of the param- 
eter a in the depth-averaged equations. At best, one may consider a as a "fit parameter" 
in the equations and the simple value a = 1 works well in our case. 

6. Results for sand 

The visual preliminary observations with sand were different than with glass beads. 
Waves were observed without forcing for slow flows but seemed to disappear for rapid 
flows in apparent contradiction with an inertial instability. In this section, we show that 
the forcing method allows to clarify the situation. 

6.1. Experimental results 

We have measured the dispersion relation of the waves for the sand with the same method 
used for the glass beads. However, important difficulties arise in the case of sand. Since 
the flow with sand is strongly unstable, it is difficult to define a linear region where 
the waves grow exponentially Moreover, the large noise associated with the natural 
instability makes the measurements much less reproducible than with glass beads. It is 
therefore more difficult to quantitatively measure the dispersion relation with the sand 
than with glass beads. In order to measure the growth rate, we have chosen the following 
method. First, the two photodiodes are located very close to the nozzle x\ — 2 cm 
and X2 = 10 cm to prevent as much as possible the non linear evolution of the waves. 
Then, for each frequencies /, the growth rate defined as a(f) = (1/(^2 — x\)) ln(a2/ai) 
is averaged over many measurements, carried out at different forcing amplitudes and 
forcing frequencies. 
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Figure 18. Results for sand 0.8 mm in diameter. Spatial growth rate as a function of frequency 
for different Froude numbers (9 = 35°). The corresponding thicknesses of the flow are h = 3.3 
mm (•), h — 4 mm (□), h — 5 mm (■) and h — 6 mm (o). 



In spite of the uncertainties and the large error bars, the general trend of the dispersion 
relation can be measured as shown in figure^] which presents the evolution of the growth 
rate cr(/) with the frequency for different Froude numbers, at a fixed angle 8 — 35°. This 
plot has to be compared with the one obtained for glass beads ffigure IT2l . We observe 
that for a given Froude number, the variation of the growth rate with the frequency is 
similar to the glass beads. The low frequencies are amplified while modes above a given 
cutoff frequency are damped. However, the order of magnitude of the growth rate for the 
sand is between 0.02 cm -1 and 0.1 cm -1 , which is about five times the typical growth rate 
measured with glass beads. This explains why the natural instability is easily observed 
with sand and not with glass beads. 

Another difference between sand and glass beads arises when studying the influence 
of the Froude number. For sand, the growth rate of the most unstable mode increases 
as the Froude number decreases by contrast with glass beads. The most unstable mode 
is observed for F — 0.25, which corresponds to the slowest flow that may be achieved 
at this angle. This result is all the more surprising since the cutoff frequency decreases 
when the Froude number decreases, as observed before with glass beads. 

These results are confirmed for other angles of inclination. Figure I19f a) gives the 
growth rate of the most unstable mode a max as a function of the Froude number for 
all the angles studied (32° < 9 < 36°). Considering the experimental difficulties, precise 
measurements are carried out for each angle only for two extreme flows: a slow flow very 
close to the flow threshold and a fast flow, corresponding to the limits of our set-up. In 
this figure, we notice that a max is much larger for low Froude flows close to the flow 
limit (grey zone) than for high Froude flows. On the other hand, a max remains positive 
even at high Froude number, which means that the flow is still unstable. Therefore, the 
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Figure 19. Results for sand 0.8mm in diameter, (o) spatial growth rate of the most unstable 
mode <T m ax{F)\ (b) neutral stability curve f c (F) for different inclination angles (given in the 
inset). The lines give the theoretical predictions (a = 1). 

preliminary observation that the waves disappear at high Froude number comes simply 
from the fact that the maximum growth rate is a decreasing function of the Froude 
number. 

It is also interesting to observe the behavior of the cutoff frequency f c with the Froude 
number (figure [T%) . We notice that f c decreases when F decreases, which is the sig- 
nature of a "high-Froude-number-instability" . However, a range of unstable frequencies 
still exists for the smallest Froude number F min w 0.25 that can be achieved in the 
experiment. This means that, for sand, the stability threshold is pushed away below the 
onset of flow. 
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6.2. Comparison with the theory 

From the measurement presented in § 4.2 we know that the velocity law for flows of sand 
has the same analytical form than the velocity law for glass beads but with different 
coefficients. The mean velocity it, the thickness h and the inclination 9 are related trough 
the following relation: 

-7 + /V- ^m- (6-1) 



y/gh h stop (0) ' 

with 7 = 0.77 and (3 = 0.65. 

From equation (|4.20|) it comes that the linear theory predicts an instability above the 
critical Froude given by: 

F.-?f 1 -S2,V (6.2) 



3 \ Vcos"# 

In the range of inclination we used, the predicted threshold is then of the order of F c 0.1. 

The stability threshold predicted by the linear stability analysis of the Saint- Venant 
equation for the sand is therefore close to zero and far below the smallest Froude number 
that is achieved in the experiment. This prediction is compatible with our experimental 
measurements (figure 119b ) showing that the flows observed with sand are unstable from 
the very onset of the flow. 

Concerning the dispersion relation, the theory does not describe the observed stabi- 
lization of the short wavelengths, as already noticed for glass beads. No cutoff frequency 
is predicted. However, we can compute for a given set of control parameters (6, F) the 
maximum spatial growth rate o ma x ■ In the theory, the maximum growth rate is achieved 
at infinite frequency. Using (|6.1|l for the law u(h, 9) and the expression of the spatial 
growth rate IjA 6J) from Appendix A, we obtain (a = 1): 



1 r / \ a 
--H+oo) = - 



(6.3) 



5o - 1 - (j) 
2(F+l) 

where a is given by (|4.18|l and Co is given by (|4.19|l . 

In figure EH a ) the lines give the prediction for (r max (F, 9). We first notice that in the 
range of Froude number where measurements are carried out, the predicted maximum 
growth rate is always positive, i.e. the flow is always unstable. Secondly, for F > 0.25, 
the maximal growth rate predicted by the theory decreases when the Froude number 
increases, as observed in the experiment. Finally, although there is no quantitative agree- 
ment between the predicted maximal growth rates and those measured experimentally, 
the order of magnitude of o~ max predicted by the theory is of the same as the order of 
magnitude as the experimental measurements. 

6.3. Discussion 

The linear stability analysis of the Saint- Venant equations seems therefore sufficient to 
understand the properties of sand flow that were first surprising in the observation of 
the unforced system. The theory predicts that the flow is always unstable and that the 
most unstable modes occur for very slow flows, close to the flow threshold. Therefore 
the occurrence of waves for slow flows in sand does not come from a new instability 
mechanism but results from the classical long-wave inertial instability. 

However, the features of the instability contrast dramatically with the instability in 
classical fluids due to the difference in rheology and flow rules. A granular media can 
be unstable from the very onset of the flow unlike classical fluid flows. This property 
actually results from the existence of a critical angle in granular flows. Unlike classical 
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Figure 20. Predicted dimensionless growth rate of the most unstable mode as a function of the 
Froude number for glass beads (6 — 29°) (a) and sand (6 = 35°) (b) (a = 1). Grey zone: no 
flow. Hatched zone: no measurements. 

fluids, the speed of the kinematic waves in a granular flow does not vanish at the flow 
threshold, i.e. when h — > h stop - If the speed of the kinematic wave cq at h — > h stopi 
Co = u(h st op) + u '\(h Btop ) hstopi is larger than the speed c + of the gravity waves, c + = 
u(h s to P ) + \J gh s top cos 6, the flow may be unstable from the very onset of the flow. This 
explains why the Saint- Venant equations can predict the formation of waves with an 
inertial mechanism even when the mean velocity goes to zero. 

Furthermore, the difference between glass beads and sand becomes clear in the light 
of the linear stability analysis. The instability mechanism is the same. However, because 
of quantitative differences in the coefficients of the flow rule, the characteristics of the 
instability differ in the range where experimental measurements are possible. This is 
clearly shown in figure 1201 where the maximal growth rate predicted by the theory is 
plotted as a function of the Froude number for both glass beads and sand. We also 
indicate the range of Froude number where experiments are carried out. It is clear from 
figure |20| that there is no qualitative difference between both systems. The difference lies 
only in the relative position of the unstable region compared to the measurement region. 
For beads the flow threshold is below the instability threshold and the measurements are 
carried out in a region where the maximum growth rate mainly increases with the Froude 
number. By contrast, for sand the flow threshold is above the instability threshold and 
the measurements are carried out in a region where the growth rate decreases with the 
Froude number. 

7. Conclusion 

We have presented in this paper an experimental study of the long surface waves 
instability for dense granular flows down rough inclined planes. By imposing a controlled 
perturbation at the entrance of the flow, we have been able to precisely measure the 
threshold and the dispersion relation of the instability. 

Using glass beads we have shown that the long-wave instability is controlled by the 
Froude number and occurs above a critical Froude number. The stability threshold and 
the velocity of the waves measured experimentally are quantitatively in good agreement 
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with the predictions of a linear stability analysis of the Saint- Venant equations. Using 
sand, we have observed that the properties of the long-wave instability are strongly 
modified: the flow is always unstable and the most unstable modes occur for slow flows, 
near the onset of the flow. Despite the apparent qualitative differences between both 
systems, we have shown that the wave formation in both cases results from the same 
instability mechanism and can be described by the stability analysis of the Saint- Venant 
equations. The difference between sand and glass beads only comes from quantitative 
differences in the coefficient of the friction law introduced in the Saint- Venant equations. 

The instability mechanism for the waves formation in granular flows therefore results, 
as for classical fluids, from the competition between inertia and gravity. However, the 
features of the "classical" long-wave instability can be modified due to the specificity of 
the friction law for granular flows. The most dramatic effect is that the stability threshold 
may be lowered below the onset of the flow, i.e. a granular flow may be unstable as soon 
as it flows, when the velocity is small. This property strongly distinguishes granular 
flows from classical fluid flows and is closely related to the existence of a critical angle 
in granular material. The existence of a strong instability close to the onset of the flow 
certainly has dramatic consequences for the non linear evolution of the waves. When 
the deformations develop, the thickness of the layer can rapidly become less than the 
minimum thickness needed to flow. The flow then evolves towards a succession of surges 
separated by material at rest. A similar behavior may be expected for other materials 
as soon as the material will present a yield stress (mud, clay, granular. . . ). This study 
then suggests that the existence of surges often observed in geophysical flows could result 
from the existence of an instability at the onset of flow due to a non zero yield stress 
condition of the material. 

Another interesting result of this study is the relative success of the depth-averaged 
equations in predicting the stability properties of granular flows. As soon as a relevant 
friction law is taken into account, quantitative properties can be predicted. This work 
then provides a good test for the relevance of the friction law deduced from the steady 
uniform flows, in a case where inertial effects determine the dynamics of the flow. 

However, our study reveals an important limitation of the depth-averaged approach 
we use for describing the instability. The simple first order Saint- Venant equations are 
unable to predict the cutoff frequency observed in the experiment. This cutoff frequency 
is the signature of a dissipative mechanism that stabilizes short wavelengths. Since dry 
granular flows do not experience surface tension, this stabilization mechanism should be 
related to the streamwise velocity variations that are second order effects in a shallow 
water description. In order to take into account these longitudinal gradients, one should 
a priori know the full three-dimensional constitutive equations, which is still an open 
problem for dense granular flows. Our measurement of the surface wave instability and of 
the cutoff frequency could therefore serve as a test for future propositions of constitutive 
equations. 
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Appendix A. 

We give in this Appendix the spatial stability analysis of the dispersion relation l|4.10|l 
given by: 



2aCbk + [{a - b)k - aCb) + [ — - a) k 



F 2 



0. 



(Al) 



The pulsation u) is real and the wavenumber k is complex: k = k r + iki, the flow is then 
unstable when kik r < 0. The resolution of the dispersion relation (|A 1J gives two spatial 
modes (+) and (— ) for k(u>) that are given by: 



k\, — 



a 



. V2(^ - ab) „ 
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where is given by: 



- 4 ( a(a - 1) + ) Co 2 



(a - bf 
F 4 



(A 2) 

L 

1 

(A3) 
(A 4) 



For the friction law considered in the article the parameter a is positive and the parameter 
b is negative. The mode (— ) is then always stable whereas the mode (+) may be stable 
or unstable depending upon the Froude number or the angle of inclination. However it 
can be shown that the sign of kik r , i.e. the stability of the flow, does not depends upon 
the pulsation Co. Therefore, to find the stability threshold we can study the asymptotic 
form of the dispersion relation. For uj — > oo, one finds: 



fc±(+oo) = 



fc±(+oo) 



a ± J a(a — 1) + 4j 



Tall 



y/ a(a - 1) + £r) 



2F* (a ± yja(a - 1) + -fr) y/a(a - 1) + J 
The mode (+) is then unstable when: 



(A 5) 
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Appendix B. 

We give here another way to determine the instability condition i|4.13[l . which under- 
lines the competition between the kinematic waves and the gravity waves in the stability 
of the flow. 

By differentiating l|4.7fl with respect to t (resp. with respect to x) and (|4.8|) with respect 
to x, the perturbed velocity field u\ may be eliminated and one obtains a partial linear 
equation for h\ given by: 

d 2 h x d 2 h x , 1 N d 2 /n a fdht , b.dhA , . 

^ + 2a ^t + [a ~^ ) ^ = -T r2 \-dt +{l --a ) ^i)^ (B1) 
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We then re-write (|B It as: 

d_ . _ d 
dt 



" + dx 



d_ 

dt 



dx 
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J2 



dhi 
~~dt 



_ dhi 
C ° dx 



(B2) 



where c± = a ± y/a(a — 1) + (1/F 2 ) is the speed of the gravity waves upstream and 
downstream and Co = 1 — (b/a) = 1 + (du/dh)o is the speed of the kinematic waves. 
Equation (|B 21) reveals waves hierarchy in the system 1 Wit ham 1974J) . Long wavelength 
perturbations propagate to the first order as kinematic waves with a velocity cq. The 
effect of higher order terms may be captured by substituting d/di ~ —dod/dx on the 
left-hand side of IB 21). which leads to: 



(go -c_)(c 4 



a 

F2 



dt 



95: 



(B3) 



This equation is a combined diffusion equation and the stability of the perturbation is 
given by the sign of the "diffusion coefficient" : (F 2 /a)(co — c_)(c + — Sq). For the friction 
law considered in this paper, a > and b < 0, i. e. Co > f . On the other hand, c_ < f 
since a > 1 and i 7, > 0. Therefore (co — c_) is always positive and the flow is unstable 
when: 

co>c"+. (B4) 
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